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Abstract 



We investigate different types of collective excitations in a quantum dot 
containing finite number of electrons at zero magnetic field. To estimate 
the excitation energies analytically we follow the energy weighted sum-rule 
approach. We consider the most general multipole excitation with angular 
momentum I, and the breathing mode excitation (monopole excitation) of a 
large quantum dot, for three different types of effective electron-electron inter- 
action. These are the logarithmic interaction, the short range pseudopotential 
and the coulomb interaction. The ground state density of the many-body sys- 
tem is calculated within Thomas-Fermi approximation. The analytical results 
for the collective excitation energies and their dependence on the system size 
and other external parameters are discussed in detail. 

PACS numbers:73.20.Dx, 72.15.Rn 
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I. INTRODUCTION 



With the advances in the nanofabrication technology, the study of the two dimensional 
electron gas(2 DEG) has become an interesting topic. Quantum dots are two-dimensional 
electrons confined within a finite area by applying various gate voltages. It is an example 
of a finite size quantum system where the number of particles can be varied from a few to 
few thousands. In recent years, these systems have received considerable attention as a tool 
to study the nature of electron-electron system in a finite size system. In particular one of 
the goals has been to determine the far-infrared response of the system and to understand 
the nature of collective excitations. Though the collective excitations have been extensively 
studied in other finite fermion systems like atomic nuclei and metal clusters |||4]], in 
three dimensions, it still remains to be analysed in detail in two-dimensional quantum dots. 

Recently, evidence has been found for a strong collective dipole mode and its splitting 
in the presence of the magnetic field in far infrared spectroscopy experiments . There is 
likely a further evidence of a quadrupole excitation in a large quantum dot in the presence 
of magnetic field |J in experiments done using far infrared spectroscopy. Recently collective 
spin and charge density excitations in a quantum dot containing 200 electrons have been 
observed @. Considerable amount of theoretical work has been done to study the nature 
of edge manetoplasmon excitations of quantum dots in strong magnetic fields @-TT]. The 



collective excitations have been studied within the classical hydrodynamical model of quan- 
tum dot HT2T[13. A lot of effort has been made to understand the nature of the collective 



excitations in these nanostructures by using different theoretical tools. 

In this paper we present detailed results on the nature of the collective excitations in 
quantum dot systems at zero magnetic field. We consider the multipole excitations and 
breathing modes of a quantum dot with finite number of electrons, and calculate analytically 
the excitation energies within sum rule approach. So far most of the microscopic calculations 
have been done numerically, and the results are therefore restricted by the particular choice of 
the device parameters. The main aim of this paper is to investigate the collective excitations 
in a dot containing large number of electrons analytically, to study the variation of the 
excitation energies with the device parameters and with the change of electron number. In 
a quasi zero dimensional nanostructure like quantum dot, the nature of electron- electron 
in teraction is not very transperent. Other than the conventional coulomb interaction, we 
consider two other types of effective interactions. The advantage of obtaining an analytical 
expression is the variation of the collective modes with the interaction strength becomes 
clear. The ground state density profile is evaluated using the simple Thomas- Fermi (TF) 
approximation. The T-F approximation gives a good account of the bulk properties of the 
system when the number of particles, N, is large \T5\ . The T-F density deviates from the full 
quantum mechanical density only in the region very close to the turning point and it neglects 
the asymptotic tail of the density. For the higher multipolarities the density fluctuations 
shift more towards the edge. But for estimating the low lying multipole excitations of a large 
quantum system, the T-F density is good enough and the edge correction of the density gives 
rise to a small correction over the leading order estimate. We estimate the edge correction 
for some important moments of density, and the corrections typically goes as, ~ c(7)/iV 3//4 , 
where the term c(7) increases with increasing multipolarity I. This estimate of the edge 
correction justifies the evaluation of different moments of a large quantum system by using 
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T-F density. So that the low-lying collective excitations of a large dot can be estimated 
within the T-F approximation. 

For the electron-electron interaction we have taken three different forms of the two-body 
interaction. These forms are chosen so as to facilitate analytical calculations of excitation 
energies. However, we also point out that in some limits these approximate the effective 
electron-electron interaction in a quantum dot. Since the effective interaction changes with 
the number of particles, it is useful to work with different approximate interactions. 

This paper is organized as follows: In section II we discuss the ground state of the 
quantum dot in the simple TF approximation assuming various forms of the electron-electron 
interaction. In Sec. Ill we discuss the RPA dynamic polarisibility of quantum dot for different 
excitation operators, and the corresponding strengh distributions of such excitations. Also 
we discuss how the collective excitation energy can be estimated using the moments of 
the corresponding strength distribution. We explicitely evaluate the multipole excitation 
energy for different models of quantum dot in Sec. IV. In Sec.V we derive the sum rules and 
excitation energy for the breathing mode of the dot, from simple scaling argument. Section 
VI contains a summary. In Appendix A and Appendix B, we outline the corrections that 
may arise from using T-F density for the ground state. 



The two dimensional Hamiltonian to describe the many electrons in a quantum dot may 
be written as 



Here we have assumed that the electrons are confined in a parabolic potential which is 
a good approximation to the device potential, m* is the effective mass of the electron in 
GaAs sample, which is 0.067 times the actual electron mass. We estimate the density 
profile of the many-body system by using T-F energy functional. It is well known that 
for the large number particles, T-F density profile agrees very well with the full quantum 
mechanical density, except very close to the turning point. Asymptotic tail of the density 
profile can be obtained by including the gradient corrections in the density functional, which 
includes the higher order h corrections in the T-F energy functional. In evaluating the energy 
weighted sum rules, T-F density gives the leading order contribution, and the small edge 
correction can give rise to the higher order (1/iV) corrections. In what follows, we assume 
that N is sufficiently large so that these corrections may be neglected. The energy functional 
describing the ground state of this many body system is 



II. THOMAS-FERMI METHOD FOR GROUND STATE 




(1) 




(2) 



where the self-consistent potential V sc is given by, 




(3) 
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Within the Thomas-Fermi approximation the kinetic-energy density is given by, 

h 2 TT 

2m* 

for spin 1/2 electrons in an unpolarised dot. The self-consistent equation for determining 
the density is obtained by minimising the free energy 

F = (E-fiN), 

where /i is the chemical potential. Thus we have the following self-consistent equation for 
the density, 



h np ^ + - m *uj 2 r 2 + / d 2 rV{\r - r'\)p{r') = /i. (4) 



- 2 _7 

m* ' 2 

The chemical potential is determined by conserving the number of particles in the dot, 

d 2 rp(r) = N. (5) 

In the following subsections we consider three different forms of the repulsive two body 
interaction and try to find out self-consistent density analytically. We first show that for 
logarithmic interaction and short-range interaction, the self-consistent density can be solved 
exactly. However, for coulomb interaction the above integral equation can not be solved 
exactly. Therefore we estimate the ground-state density profile by variational method. 



A. Logarithmic interaction 

The two dimensional Hamiltonian with logarithmic interaction is given by, 

h = h£r. + -ffEM^l- (6) 

i=i z " i z i<j u 

Here, a is an arbitary parameter which determines the range of the two-body interaction, 
and g is the strength of the interaction. The advantange of taking logarithmic two-body 
interaction is that the Thomas- Fermi model is exactly solvable JTjyiTj, and it can mimic the 
realistic two body interaction. This model may also be justified from the finite thickness 
effect of the two dimensional system |17|]. As shown by Zhang and DasSarma ||18||, the 



effective interaction for a quasi two-dimensional electron system in momentum space is 
given by, 

e 2 1 + 9k/8b z + 3k 2 /8b 2 z 

where, b z has dimensions of inverse length. The effective potential given above behaves 
approximately as log{r) over small distances and - for large distances. The strength of the 
logarithmic interaction can be written as g ~ ^p*. 
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Following Thomas-Fermi method the density distribution of the electron system can be 
determined from the following equation, 



TTfl 2 p(r) 1 * 9 2 f 72 / A 7 \^~ T '\ 

+ —m to r —gar p(r )ln\ 1 = p. (8) 



m* 2 

Now applying laplacian operator on both sides of the equation we obtain the following 
differential equation, 

—V 2 p(r) + 2m* J 1 - 2ixgp{r) = 0. (9) 

Assuming the circular symmetry for the density of the ground state, the above equation can 
be writen as, 

d 2 y Idy 

-7-2 + -~r ~ 2gy = °' 10 ) 

ax z x ax 

where, the dimensionless variable x = r/l and y = p(r) — ^p^. The harmonic oscillator 



length scale and the dimensionless coupling are defined as, l Q = and g = The 

solution of the above differential equation is given in terms of modified Bessel functions, 



y(x) 



AI (j2 Z gx)+BK (j2g~x) J (11) 



where A and B are constants to be fixed by the boundary conditions. In the limit x — ► 0, 
Jo — > 1, where as K — > —ln(x) [|1^]. The latter condition implies that the density at the 
centre is infinite which is unphysical since it would require infinite energy to push particles to 
the centre due to the Coulomb repulsion. Obviously minimisation of the ground state energy 
implies that B should be zero. The quantity A is then fixed by the boundary condition, 
p( r o) = at the turning point. We have therefore the full solution for the TF density, 

p{r) = -—{1 - , forx<x . 12 

TiVq MJ2qxn) 



7rl 2 g I (y/2gx ) 
The turning point r can be determined from the normalisation condition, 

N = f p{r)d 2 r 



1 { L 2 L h(L) 

~9 2 ^ hiiy 



, - L ^} (13) 



here L = ^/2gxQ } and x = r/r . For L >> 1, (which is equivalent to large N limit) the 
density becomes almost constant at the value p = an d the relation between L and the 
number of particles N becomes, 

2gN = L 2 - 2L. (14) 



In the other limit L << 1 (which is equivalent to the weak coupling limit), the density 
becomes almost parabolic, and the dependence of L on N is given by, 
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y 16 V 12 ; 

N=\xt{l+ l -gxl). (15) 

In fig. 1 we have compared the scaled density profile of the electrons with logarithmic inter- 
action with the parabolic one, for two different values of the parameter L. When L — 1, 
then the density becomes almost parabolic, but for L = 20, the density becomes very flat. 



B. Short range interaction 

We now consider the case when electrons in the dot are interacting by a very short range 
delta function interaction, V(r) = VoS(r). For a dilute Fermi gas the two body interactions 
can be approximated by delta function pseudopotential. The interaction strength Vq is 
related to the s-wave scattering length a in the following way, Vo = 4 ^"» a ■ Also it has been 
found that for 2-D electron gas in strong magnetic field the short range pseudo-potential is a 
very good approximation . This system is also interesting because short range interaction 



among the fermions can mimic exclusion statistics f2T|j . 

The energy functional of the electrons with short range two body interaction can be 
written as, 



E[p) = J d 2 r[^P 2 (r) + ^r) + ^p 2 (r)]. (16) 

By minimizing the energy functional with respect to the density, we obtain the following 
equation, 

h 2 7i 1 

—p(r) + -m*uj 2 r 2 + V p(r) = /i, (17) 
m* 2 

where \x is the chemical potential of the system. The self consistent density can be written 

as, 

P( r ) = ^(^-^ 2 )> ( 18 ) 

where, x is the dimensionless variable x = t/Iq and other parameters are defined as, a s = 
1 + V °™i and ft = -j^. The turning point of the density distribution can be determined from 
the conservation of the total number of particles, 

N — [ d 2 rp(r) = (19) 

The results thus obtained with short range interaction are similar to the non-interacting 
but confined electrons. The only difference with the non-interacting electrons is the scaling 
factor l/a s both in density and in the number of particles. 
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C. Coulomb interaction 



Next we consider the case where the two-body electron electron interaction is 1/r. For 
this form of the two-body interaction, T-F self consistency equation for density can not be 
solved analytically. To do the rest of the calculations analytically, we take a variational 
ansatz for the self consistent density profile and determine the parameters of the ansatz by 
minimising the energy functional. In zero magnetic field we can take the following ansatz 
for the density profile, 

^) = ^4-(zo-z 2 ), (20) 

where, a c is the variational parameter, and x is the dimensionless variable r/l . The density 
of the non interacting electron gas corresponds to the limit a c — 1. The inclusion of the 
two body repulsion flattens the density profile, which can be thought of as the reduction of 
the frequency u> e ff of the effective mean field potential. So the electron-electron interaction 
makes a c > 1. In the case of short-range interaction the density profile is parabolic, and also 
for the logarithmic interaction density profile becomes parabolic in the strong confinement 
limit. We expect at zero magnetic field the ansatz can give good estimate of the ground state 
density in the strong confinement limit. We must, however, caution that this ansatz breaks 
down in the presence of a strong magnetic field. The ground state density here deviates 
from the parabolic form and is approximately p(r) ~ J 1 — (r/R) 2 p2] . For convenience 
we define p = 27r /s Q ■ The turning point x can be calculated from the total number of 
electrons N, 

N = J d 2 rp(r) = ^ P o4 (21) 

Let us calculate the contributions to the energy coming from the different parts of the density 
functional as a function of the parameter p. Contribution to the kinetic energy is given by, 

E T = — [ d 2 rp 2 (r) = ^ 7r2p ° / ° 2x ° . (22) 
2m* J ' 6m* 

The potential energy due to the external confinement is given by, 

E p = ^m*u 2 J d 2 rr 2 p(r) = -^7rm*u 2 p l*x 6 . (23) 
The interaction energy within the Hartree approximation is given by, 

E H = ^J d 2 rd 2 r'p{r)V{\f- r'\)p(r), (24) 
where, V(\f — r'|) = ^ A ■ In momentum space the above expression can be written as, 

E » =1 2&eI d2 V^-^ < 25) 

where, 
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p(q) = f p(r)e^d\ = 4w £^. (26) 
J q 

After performing the integration we obtain the Hartree energy in the following form, 

r(|)r(4)rg e 2 . 2 2 

Eh ~ 2r(5/2)r(9/2)r(5/2)7 (7rx ° Po) ■ (27) 

In terms of the variational parameter a c and the number of particles N, the total energy of 
the system is given by, 

Eo = %oNn*l /2 + 4*) + l^^l^^ 4 - ( 28 ) 
<J etc -J-L^Trac' e< o 

To get the minimum energy, the total energy is be minimized with respect to the variational 
parameter a c when the number of particles N is held fixed. The a c that one obtains at the 
minimum of the energy is therefore a function of N, among other things. This is in contrast 
to the short-range interaction where a s is independent of N. 

Having thus determined the ground state energy and the density in the TF approxima- 
tion, we now use these in evaluating the collective excitation energies for different types of 
collective excitations. 



III. SUM RULES AND COLLECTIVE EXCITATIONS 

The excitation spectrum of any system is usually probed by applying external fields. The 
effect of the perturbing probe can be described through an interaction Hamiltonian, 

H int = X(Fe iut + F^e- iut ), (29) 

where the quantity A gives the strength of the oscillating field. Given an excitation operator 
F, many useful quantities of the perturbed system can be calculated from so-called strength 
function, 

S± = E I < n \ F ±\® > \ 2 &{E - E n ), (30) 

n 

where, E n and \n > are the excitation energy and excited state respectively, and F + = F, 
F_ = Ft. Various sum rules are derived through the moments of the strength function, 
defined as, 

m $ = \j ' E\S + {E)±S_{E))dE 

= i(< 0\F{H - E ) k F^\0 >±< 0\F\H - E ) k F\0 >). (31) 

It is easy to see that, for a given k, the moments may be expressed in terms of the com- 
mutators of the excitation operator F with the many body Hamiltonian H . We give below 
some of the useful sum rules, 
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m, 



,! 2 
1 



-<0|[Ft,F]|0>, (32) 



m 



m 



2 



<0\[F\[H,F])\0> (33) 
<0|[J f ,J]|0> (34) 



mj = ^<0\[J\[H 1 J}\0> ; J=[H,F}. (35) 

The choice of the excitation operator F is dictated by the physics that one wants to 
describe. In this paper we consider two types of collective excitations: the multipole exci- 
tation modes and the breathing modes in a the quantum dot. In two-dimensions multipole 
excitation operator can be written as F = Y^i r\e * = J2i z \i where z is the complex variable 
x + iy. Angular part of any function in two dimensions can be expanded in terms of these op- 
erators. These multipole modes describe the angular momentum excitations of the system. 
The excitation operator for breathing mode is given by F = J2i r 1- The collective excitation 
energy can be estimated from the sum-rules of the corresponding excitation operators. 

For a highly collective state the strength distribution becomes a sharply peaked function 
around the collective excitation energy. If the excitation operator F is hermitian, or if the 
Hamiltonian does not contain any time-reversal symmetry breaking term (as in the quantum 
dot without magnetic field), then we can take = and the low lying collective excitation 
energy is given by, 



E, 



\ 



c 

mi 



(36) 



Near the collective excitation energy, we can approximate the strength distribution by a 
delta function S±(E) = a±5(E — E c ), then we can trivially show the above form of collective 
excitation. But in general the strength function has a finite width. Following ref. 0, one 
can then derive the above form of collective excitation energy for multipole excitations by 
using the variational principle. Given the N electron ground state |0 >, it is possible to 
find the collective excitation energy and the collective state |c >, if one is able to find an 
operator 0\ which satisfies the following equation of motion, 

[H,0^}=huj co uOl (37) 

The state O^\0 > has excitation energy hu co u. The excitation energy is then given by the 
following expression, 

<0|[O,[g,Ot]]| 0> 

huJcM - <o|[o,ot]| > • (38) 

We may now take the variational ansatz for O* as, = F + aJ with a as the variational 
parameter. Substituting this expression in the above equation, we obtain the excitation 
energy as, 



f^coii = — 39 
zamj 
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where we have assumed m 2 = 0, for the system in the absence of the magnetic field. By 
minimizing the energy with respect to the variational parmeter a, we obtain the collective 

excitation energy as hw co u = \ = E c . 



IV. MULTIPOLE EXCITATIONS IN QUANTUM DOT 

The multipole excitations in a quantum dot lead to the change in total angular momen- 
tum. In particular, the dipole excitation in parabolic quantum dot has been well studied 
both theoretically and experimentally. The corresponding operator for the dipole excitation 
operator is given by F = x + iy. One can generalise this and write for any multipole exci- 
tation operator Fi = z l . We use this form to calculate the relevant sum-rules for multipole 
excitations. The first moment of the strength distribution function may be written as, 

rof(0 = ^<oiE4[£,E*> 

~J dW'Mr), (40) 

where p(r) is the ground state density distribution. 

The calculation of m^(l) is more complicated. To make the analysis transparent, we 
consider the various contributions to the third moment separately. The sum-rule m^(l) may 
then be written in the following form, 

m$(l) = m 3 (T) +m 3 (V) + m 3 (ee), (41) 

where the first term denotes the contribution from the kinetic energy operator in H, the 
second term denotes the contribution from the confinement potential, and the last term is 
the contribution from the interaction part. We now calculate separately each one of these 
contributions. 

The kinetic energy term in m 3 (T) is given by, 

m 3 (T) = 2h\l - / d 2 r[r 2(l - 2 \lr(r) + 2(1 - 2)A(r)], (42) 

where, r(r) =< >, and A(r) =< J2i ^172 >=< % m *}l >■ Here the angular brackets 
denotes the averages over the momentum distribution. Note that unlike the ground state 
where only the first term contributed, for the multipole excitations the centrifugal barrier, 
denoted by the second term, contributes to the sum rules. Using the Thomas- Fermi approxi- 
mation this contribution is givenby, A(r) = \r(r) = f^rP 2 ( r )- Substituting these expressions 
in the above equation we obtain, 

ms(T) = 2 " 4 ^; 1 ) 2 g / dV-VM. (43) 

The contribution coming from the harmonic oscillator cofinement in the sum-rule is given 

by, 
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m 3 (V) = lm*ulh 2 ^ J d 2 rr 2 ^p(r) 
The last term 777.3 (ee) due to two body interaction can be written as, 

m3 ( ee ) = -h 2? ^[ f d 2 rV H (r)r 2 ^ p" (r) + (2/ - 1) f d 2 rV H (r)r 2l - 3 p'(r) 

+ Jd 2 rJ d 2 rV(r)r'- 1 e- a V(|f-r'|)e iW V"- 1 p'(0] ) (45) 

where V(r) is the two body interaction, Vh = J d 2 rV(\f — r'\)p(r'), and p'{r) = ^r^. 

We thus have the general expressions for the sum rules where the only unknown is the 
ground state density p{r). For specific interactions, the ground state density was calculated 
in Section II. We now use this to explicitly compute the sum rules. Note however, that 
the T-F density calculated in section II fails near the turning point, since the T-F density 
abruptly goes to zero. But for the lower multipolarities main contribution to the moments 
comes from the bulk region, and the diffusive edge rigion gives rise to small correcction, 
which decreases with the increasing N. We discuss the corrections to this approximation in 
the appendices. 



(44) 



A. Logarithmic potential 

For logarithmic potential the first moment mf is given by, 

. ,h 2 i 2 



777; = 777 



"1 ~~ *2 

m 

.2/2 



(frr 2 ^ p(r) 



= m'^pjfr* fdxx 2 ^[l - ^H], (46) 

777*^ Jo I (L) 

where the central density p c = — ^| and the dimensionless variable x = r/r , and L = ^j2gxQ. 

In the case of the third moment 7773", kinetic energy and the confinement energy contri- 
butions can be written as, 

777 3 (T) + 777 3 (V) = ?%[2(i - I) 2 — [ d 2 rr 2 ^p 2 (r) + lm*u 2 h 2 [ d 2 rr 2 ^p(r)]. (47) 



777 777 



We now evaluate the remaining term 7773 (ee). Integrating eq.([|5|) by parts we can write 
7773 (ee) as, 



til r r r)V 

m3 ( ee ) = — — [ / d 2 rV 2 V H r 2 ^p(r) + 2(1 - 1) / d 2 r^-r^p(r) + J 3 ], (48) 
2m* z J J or 

where, Vh(t) is the hartree potential. For logarithmic interaction we can use the following 
relations, 

Wtf(r) = -2irgp{r), (49) 
9V H * 2 h 2 ndp 

— — = -777 u r - — . 50 

or 777* or 
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Therefore the term ^3 is given by, 

J 3 = J d 2 r J d 2 rp(r)' r l ~ x e~ lW V{\r - r'\)d w r' l ~ x p\r') 

= (2tt) 2 f d 2 qf{q) f drr 1 J far) pf '(r) f dr'r' 1 : Ji{qr')p' '(/), (51) 



where /(g) is the two body potential in momentum space, and p'(r) = After doing 

the integration by parts and using some identities for Bessel functions we can rewrite the 
above expression in the following form, 



d 2 qq 2 f{q) I d 2 re 1 ^ e~ l ^ e r {1 ^ p[r) J d 2 r' e^' e l ^ e ' r'^ p(r') 
= 2ng J d 2 rr 2 ^p 2 (r). (52) 

Using the above results we can write the multipole excitation energy in the following form, 

„ mt - , 3ft 4 7r(Z - l) 2 Jd 2 rr 2 ^p 2 {r) 
mj m A J d A rr^ l ~ l > p[r) 

The above equation can be written as, 

E 2 c =hW[l + 6l(l-l)^B], (54) 

where, 

l_l Jo 1 dxx 2l - 3 [l - 

~ ] ii^x^[i Io{L) 

In the large L limit, one can approximate the function [1 — by a simple form [1 — 

exp(L(x — 1))], and the approximate value of f(L) is given by, 

~ = 1 - 2F 1 (2l -2,21- 1; L)e~ l + F l {2l - 2, 21 - 1; 2L)e' 2Z 
1 J l-F 1 (2Z,2Z + l;L)e- £ ' 1 J 



^ 'V T.^-.n^r (55) 



where F\ is the degenerate hypergeometric series ||19|| . In the large L limit asymptotic N 
dependence of the collective excitations are given by, 



E c = Wl + ^. (57) 



3l(J_ 
T 2 N 

By substituting 1 = 1, for dipole excitation we can see that the excitation energy is hu, 
which obeys Kohn's theorem [23]. In the other limit L « 1, (which is the weak coupling 
limit) the dispersion relation of the multipole modes are given by, 

E c = hcu[3l - 2 - Q^j-^-~gN 1 l 2 } 1 / 2 . (58) 

We can estimate the strength of the coupling from the relation g = given in sec. II. In 
GaAs sample the values of the paraters are b^ 1 = 58A and e = 12.6. For tkv = 5.4meV 
our estimate of the dimensionless coupling is g ~ 1.36. For this value of the coupling the 
variation of excitation energy with number of electrons for 1 = 1 — 4 are shown in figure 2. 
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B. Short range interaction 



Now we calculate the sum-rules for the multipole excitations in the case of the short 
range interaction. For short-range interaction the density profile is given by, 

«0 ft s 2 

where, a s = 1 + ^p-, and /q = The first moment mf can be written as, 

h 2 i 2 



m+ = m* — - I d 2 rr 2(l 1} p( r ) 
m* 2 J 



h 2i2 ;2(i-l) 2(1+1) 

= m ^ ( 2^(/ + l) ) - (60) 

Now let us calculate the third energy weighted moment mf. Different contributions of this 
sum-rule are given bellow, 

= 4 /, -i\2*2,2 



m 3 (T) = / d 2 rr 2 ( z ~V(r) 

771 771 J 



h\{l-l) H 2 l 2 2( ;_ 3 ) 2(1+1) / fil x 

m*«2/(/ + 1) m* 2 1 J 

,fr 2 Z 2 



m3 (y) = l m *hW^ [ d 2 rr 2 ^p(r) 
m* J 



h 2r2 ,2(1-1) 2(1+1) 

= i m *ftVg- V° *° - (62) 
m* 2 2a s /(/ + 1) v ; 

m 3 (ee) = 0. (63) 
From the above expressions collective excitation energy is given by, 

E 2 c = ^r = IhW + 2(1 - 1)— . (64) 
m{ a s 

Once again for 1 = 1, the excitation energy is E c = hou consistent with Kohn's theorem. 
For the value of the dimensionless parameter a s = 3.5 the dispersion relation of the modes 
/ = 1 — 4 are shown in figure 3. 



C. Coulomb interaction 

For coulomb interaction we have taken a variational density of the following form, 

P(r) = po{x 2 -x 2 ), (65) 

where po = 2 nt 2 a anc ^ x = r /^- Here a c is the variational parameter determined by 
minimizing the ground state energy as in Section IIC. This density is similar to that of short 
range interaction. The first moment therefore mf is given by, 
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fj2 /2 __ 2(m) 

+ _ n 1 Wo (m) 

The sum of first two terms of the third moment m%(T) and 777,3 (V) can be written as, 
m 3 (T) + m 3 (V) = — [ ltm%l + l) + m - 1} ]■ (67) 



We now calculate V$(ee) for 1/r interaction, which can be written as, 

2^S [ / ^VV^-)p(r) + 2(Z - 1) / rf 2 r^ 



m 3 (ee) 



+ d 2 r dVr<-y(r)e-^^^e l 'V<-y(r')]. (68) 
J J \j? — 7™' 

where // = The hartree potential V#(r) is given by, 
y^(r) = - / dy_I_ p ( r ) 

6 J — 7*' J 

^<> 3 iz§ F(1/2 '- 3/2 ' 1;( ^ ) ' (69) 

After performing other integrals and some algebra we arrive at the final expression for 

777,3 (ee), 

The collective excitation energy for the multipole modes are given by, 

Of, 2 , i 2 

£a( i ) = [( i _i)f^L + /nV 

e 2 {AN) 1 ' 4 

E c (l) = fuv[2{l - l)/a c + 1 + 1(1 + x 

[g£|g - ^ ,^(3/2, -1/2, + 2,2:1)]]*, (71) 



where 3^2(01, a2, a3; 61, 62; x) denotes the generalised hypergeometric function [19] . The 
validity of Kohn's theorem is easily checked. For the dipole mode I = 1, we have 
3 F 2 (3/2,-l/2,2;3,2;l) = F(3/2, -1/2, 3; 1) = r{ lf$$ /2) [0- Substituting this into the 
above equation we obtain the dipole energy, E c (l) = fwj. 

In ref. |§ the multipole modes of a quantum dot with hu = 5.6meV containing 56 elec- 
trons has been calculated by using self consistent Hartree-Fock calculation. The multipole 
excitation energies for I = 2 — 4 are given by, 7.5, 9.0, 10.4 meV respectively. From our 
analytical expression we find the corresponding multipole excitation energies are 7.5, 8.9, 
10.12 meV, which are in good agreement with the self-consistent numerical calculation. For 
hu) = 5AmeV the multipole modes for 1 = 1 — 4 are plotted in the figure 4. Except in 
the case of dipole mode, the collective excitation energy decreases with increasing particle 
number. 
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V. BREATHING MODE EXCITATION 



In this section we calculate the sum-rules and the excitation energy of the breathing 
mode in a quantum dot. This excitation is also known as the monopole excitation. In the 
breathing mode excitation average radius oscillates around its unperturbed value with a 
frequency of the excitation energy 0. In a recent experiment on Bose-Einstein condensate 
this phenomena has been obseved very clearly [p4|. The monopole excitation operator is 
given by, 

F = J2^ + y")- (72) 

i 

The first sum-rule mf can be written as, 

mi = ^<0\[F\[H,F]]\0> 

= ^<0|£r?|0> (73) 
m* ^ 

which is nothing but the mean-squared radius in the ground state. 

The third moment m 3 of the strength distribution can be calculated by using a sim- 
ple scaling transformation of the ground state wave- function |p5|| . Consider the following 
transformation of the ground state wave function, 

|0 r? >=e r ' J |O>, (74) 

where J is an anti-harmitian operator J = [H, F]. Now J can be written as, J = -^-(1 + 
x If we take a new set of coordinates x' = log(x) and y' = log{y) 1 then the operator 

J becomes + -£j + Using these relations, it can be shown that the wave function 

and the coordinates undergo following scaling transformation, 

<j>r,(x, y) = e r '(j) (e v x, e v y), (75) 

where, fj = —^-77. The sum-rule 7713 can be calculated from the scale transformed wave 
function (j) v , in the following way, 

m 3 = i<0|[A[#, J}}\0> 
1 cP' 

= W < ^ > |r?= °- (76) 
Under this transformation density scales as, 

p,(x,y) = e%(e^,eV (77) 

We calculate the sum-rule m 3 seperately for different terms of the Hamiltonian. The sum- 
rule m 3 can be written as, 

m 3 = m 3 (T) + m 3 (V) + m 3 (ee). (78) 
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First let us calculate the contribution m 3 (T) coming from the kinetic-energy density. Within 
the Thomas- Fermi approximation we have, T[p] ~ / p 2 {r). So the kinetic energy of the scaled 
density is given by, 

T[ Pv ] = <?*T\p\. (79) 
Using above relations m 3 (T) can be evaluated in the following way, 

= ^?n P ]. (so) 

Similarly, for qudratic confinement potential, the potential energy scales as V v = e~ 2f, V ', and 
m 3 (V) is given by, 

m 3 (V) = 2 A V. (81) 

lit 

Finally we can calculate the term m 3 (ee) due to the electron-electron interaction, in the 
following way, 

m 3 (ee) = \^\ J d 2 r J d 2 r'p(r)V( l -^^)p(r% (82) 

where a = e**. We now calculate the above quantities and monopole excitation energy for 
different types of interaction. 



A. Logarithmic interaction 

Using the density profile of the particles interacting by logarithmic potential, the first 
moment m\ for monopole excitation can calculated exactly, 

mi = / d 2 rr 2 p{r) 

m* J 

m* 1 I (L) K L L 2 n 

5o = 

the third moment m 3 is given by, 



where the central density is po = and L = y/2gx . The sum of the first two parts of 



;2 fc 2 



m 3 (T) + m 3 (V) = 2(— ) 2 [^ J d 2 rp 2 + -m*co 2 J d 2 rr 2 p(r)]. (84) 
The third term m 3 (ee) is given by, 

1 9 If (' \T T\ 

m ^ ee ^ = 2drf^~2J J ^ r ' p ^ log ^ & a M r )^=° 



drf 

[log{a)N 2 /2] = 0, (85) 



i d 2 



2 drf 
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where a = e 2h2ri / m * . The monopole excitation energy is given by, 



i& = 2ftV[l + ^S], (86) 



where f m (L) is given by, 



UL) = \ * [ ~ , (87) 



2 /o ife 2[1 _M^) 



J 

Using the approximate form for the density profile for L » 1, the function f m (L) is given 
by, 

f rf x l-2F 1 (2,3;L)e- Z + F 1 (2,3;2L)e- 2£ 
1 - F 1 (4,5;L)e- i 
The asymptotic N dependence of the breathing mode is given by, 



E m = huj ] j[2 + ^}. (89) 

In he weak coupling limit (L « 1), the breathing mode excitation energy is given by, 

E m = 2hu>[l-lgN 1 ' 2 ] 1 ' 2 . (90) 

o 

For g = 1.36, the breathing mode is shown in figure 2. and the excitation energy becomes 
almost constant at the value \[2Tuo. This mode has higher excitation energy than the other 
multipole modes. 



B. Short range interaction 

The form of the short range pseudopotential is V = Vo5(r), and the first moment mi for 
the breathing mode is given by, 

2h 2 



mi = / d 2 rr 2 p(r) 

m* J 



6m*a 

where a s = 1 + ^pr, and x = r/l . After doing some algebra the third moment m 3 can be 
written as, 



fi 2 l 2 

° xg, (91) 



m 3 = m 3 (T) + m 3 (V) + m 3 (ee) 



m* 7 L 2m* v 12vr^a/ 2 v 12a/ 2 v 127r/ga/ J 
m* ; v 12a s ; v ; 
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The excitation energy of the breathing mode is given by, 



E=J^ = 2hu. 
V nil 



(93) 



Therefore for the short range interaction the excitation energy is independent of the inter- 
action strength. 



VI. COULOMB INTERACTION 

Now let us consider the case where the electron-electron interaction is 1/r coulomb 
interaction. Variational density of the N electron system is taken to be parabolic p = 
Po( x o — x 2 ). Where the central density is p — 2n ^a • ^ e first sum-rule m± is given by 

,2h 2 



mi = ( ) [ d 2 rr 2 p(r) 

m* J 

nh 2 Polfat 



3m* 



(94) 



Using the scale transformation properties of the density we obtain the second sum rule 7713 
in the following form, 

1 2h 2 

m 3 = -(—) 2 [A(T[p] + E p ) + E H ] 

Zi ill 

l,2fi\ 2r 2 5 ,rfft 1 . r(4)r(l/2) ,e 2 l 2 , 4 2 7l 



From the above expressions of the sum rules, the monopole excitation energy(y m 3 /m!)can 
be written as, 

F 2 - 2^ 2 rl 4- 1 4- 3r(4)r(l/2) e 2 x 
[ a c + + A(T(5/2)) 2 T(9/2) [ d hu; ) aJ , 

In Fig. 4 the N dependence of the breathing mode is shown for huj = bAmeV and it remains 
constant at the value 1.7hu. 



VII. SUMMARY 

In this paper we considered the linear response of quantum dot with finite number of 
electrons under the time dependent perturbation. We have considered different types of 
collective excitations of the dot at zero magnetic field. For a suitable excitation operator F, 
the imaginary part of the dynamic polarizability (or strength distribution) which is given 

by, 
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S±(E) = J2\ < 0\F±\n > \ 2 5{E-huj n ) (97) 

n 

becomes a highly peaked function near the collective excitation energy. The collective ex- 
citation energy is well estimated by the energy weighted sum-rule approach, and is given 
by, E co u = ^J~^- We have disscussed in general both multipole excitations(also known as 

edge modes) and breathing mode excitations in a quantum dot. To calculate the collec- 
tive excitation energies we have taken three different types of electron-electron interaction, 
which are, the logarithmic interaction, the short range pseudopotential and the coulomb 
interaction. Within the Thomas-Fermi approximation the ground state density has been 
calculated. This is a good approximation for large number of particles. For multipole ex- 
citations the angular momentum of the system changes, and the the excitation operator 
is defined as z l , which is the basis for the expansion of the angular part of any density 
fluctuation in 2-dimension. In quantum dot with large number of electrons and logarithmic 
two body interaction, multipole excitation energies form a band like structure around the 
energy hu, and the energy spacing goes as, ~ l/ r o- For short range interaction excitation 



energies are E = hujyl + 2(1 — l)/a s , which is independent of the number of electrons N, 
but depends on the interaction parameter a s . For coulomb interaction excitation energies 
take the following form, 

E 2 = nW[l + (l-l)/a c 

+ «< + ^i^W^) - (TTT)^ (3/2 ' ^ < + 1; < + 2 - 2; ^ < 98 > 

For dipole mode all results are consistent with Kohn's theorem, and the dipole excitation 
energy is Hcu. 

For the breathing mode excitation, average value of the collective coordinate r 2 oscillates 
around its equilibrium value. For logarithmic interaction breathing mode excitation energy 
becomes \f2huj in large N limit. For short range interaction excitation energy is 2hu, which 
same as the excitation energy of the non-interacting system. For coulomb interaction the 
breathing mode remains almost constant at the value l.lhu. As shown in the Appendices, 
the low lying excitations are well described by the T-F approximation for the ground state 
density in the large N limit. The corrections may also be calculated in order to correct for 
the inaccqurecies that may arise for higher excitations. 

Apart from the single particle excitations, collective excitations in large quantum dots 
and their dependence on the number of particles and external magnetic field has become 
very interesting topic. Detailed study of the collective excitations can give information about 
the electron-electron correlation in finite size quantum dot. 

I would like to thank M. V. N. Murthy for his helpful comments and critical reading of 
the manuscript. I would also like to thank R. K. Bhaduri for his helpful comments. 



APPENDIX A: 

In this appendix we describe the corrections of the collective frequencies due to the 
edge corrections of the T-F density for a simple case. Since for the logarithmic two-body 
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interaction, the T-F density can be written in a closed form, we take this example. The 
density distribution for logarithmic interaction is given by, 

p(r) = -L[l- J °/f^U (Al) 

^ ; irl 2 g [ I (V2Tx r 1 ' 



where, x = r/lo, and ro is the turning point. When L = y/2gxo >> 1, then the den- 
sity becomes almost flat. In this limit we can approximate the T-F density by a "Fermi- 
distribution" like density, which has a diffusive tail also and this density is given by, 



Po 

1 + exp(^) 



Pair) = ; , .Zr-RV (^2) 



where, p = R ~ r , and a ~ l /y/2g. The multipole excitation energies are given by, 

2 _ 2 2 3fiSr(l-l) a Q a (Q 

Ec " ^° + — ^guo — ' ( } 

where, (^2(0 — / r 2l ~ 3 p a (r) 2 dr, and Qi(Z) = / r 2l ~ l p a {r)dr. Now we can expand Qi(Z) and 
Q 2 (0 in a series of a/R, which is known as leptodermous expansion. Upto the leading order 
correction Q\ and Q 2 are given by 



Q 1 (l) = ^[l + 21(21 -1)^)% (A4) 

Q 2 (0 = ^^[l-2(/-l)(|)], (A5) 

if we set a = we recover the results given in section IV. Using the above results, upto the 
leading order in a/R, collective frequencies are given by, 



r, + r / 3/(/-l) 6/(/-l) 2 n 



So the leading order correction goes as, ~ (//viV) 3 , which is negligible for low lying excita- 
tions in a large system. 



APPENDIX B: 

In this appendix we estimate the correction of two important moments of the density, 
due to the diffusive edge of the density. These moments essentially controll the excitation 
energies we are interested in. To include the asymptotic tail of the density, we write the 
total density in the following way, 

p(r) = p TF (r) + 5p(r), (Bl) 

where Ptf{ t ) is the Thomas- Fermi density, and 5p is the edge correction. Since the quantum 
mechanical density of the electrons in the harmonic oscillator potential asymptotically decays 
as ~ e~ r ' o, we choose the following form of the edge correction, 
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5p(r) = 5p e( r °- a ) 2 /V r2/i ° for r > r - a, (B2) 

where, a is the size of the region near the turning point ro where the T-F density deviates 
from the actual quantum mechanical density. In the large N limit a « r .At r = r — a 
we match two densities, Pt-f(J'o — a ) = 8po- Now we consider the following moments, 

h(n) = J r n p(r)d 2 r (B3) 
J 2 (n) = J r n p 2 (r)d 2 r. (B4) 

The corrections of the moments due the asymptotic tail of the density distribution may be 
written as, 

poo 

Shin) « 5p e r "> /l o / d 2 rr n e- r2/l ° (B5) 

Jro 

5I 2 (n)^(5p ) 2 e 2r o /l '> f°° d 2 rr n e - 2r o /l o . (B6) 

Jro 

By using parabolic T-F density, we can calculate (h(n)) T F and (h{n))TF- The leading order 
fractional change of the moments due to the edge correction can be written as, 

6h(n) . (n + 2)(n + 4) a. 1 



=3 



(h(n)) T F 2 l r 

_ (n + 2)(n + 4) a 1 

2 V (4iV) 3 / 4 ' 1 J 

where ro = tq/Iq, for noninteracting electrons ro = {AN) 1 / 4 . The correction of h(n) is given 
by, 

51*2 (n) , , , . , . . , N/ a N r,l r/ a N 1, 
' (n/2 + l)(n/2 + 2)(n/2 + 3)(-) 2 I 3 [(4-) + H 



{h{n)) T F ' ~ ~ ~ k ''n' Al '0 

- <"/ 2 ) 3 2 [(4 !> + ^ (B8) 

Therefore the errors are at most of the order of n 3 /iV 3//4 . We may now estimate a/lo in 
the following way. The wave length A of the particles near the fermi surface is of the 



order of ~ p F ( r ) - ^he l° ca l fermi momentum is given by, Pp(r) = y2m(/i — V(r)). The 
semiclassical approximation breaks down when ~ 1. From this condition we obtain 

a /h ~ —tta = 1/(2(4A^) 1 / 12 ). This estimate shows that the error can be reduced to even 

2r 

smaller value in the large N limit. 

A more direct way of estimating corrections to the T-F limit is to compare the results 
with those obtained using the quantum mechanical density. This is rather hard in the 
case of interacting system, but easier in the non- interacting limit. The exact quantum 
mechanical density can be written as a sum over the squares of various orbital wavef unctions, 
p( r ) — En,i |Vv( r )| 2 5 where the sum is restricted so that the number of orbitals included 
conserve number of particles N. We calculated the excitation frequencies for two values of 
N, by using quantum mechanical density, and compared them with frequencies obtained 
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from T-F density. For N — 12 and I = 4 excitation frequency T-F model deviates by 
only 10 percent, for the lower frequencies error is much less. For N = 20 and / = 4 
collective mode the error due to semiclassical approximation reduces to 6 percent. This 
simple analytical calculation supports our large N asymptotic analysis of the error due 
to the use of semiclassical density. If we consider the low-lying collective excitations in a 
large quantum system, then our calculation indicates that the error in the collective modes 
decreases asymptotically with the increasing number of particles. 
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FIGURES 



FIG. 1. Comparision of the density profile of electrons with logarithmic interaction with the 
parapolic density profile for two different values of L. Solid line shows the density of the electrons 
for L = 1, dashed line shows the density for L = 20, and dot-dashed line shows the parabolic 
density. 

FIG. 2. Excitation energies (hu c ) as a function of the number of particles (N), for logarithmic 

interaction with g = 1.36. Solid line shows the exact result and the dashed line shows asymptotic 

value of the excitation energy. I values indicate the corresponding multipole modes and n = 1 level 

indicate the breathing mode. 

FIG. 3. Excitation energies (huj c ) as a function of the number of particles (N), for short range 

interaction with a s = 3.5. Multipole modes are denoted by corresponding I values and n = 1 level 

indicate the breathing mode 

FIG. 4. Excitation energies (hu c ) as a function of the number of particles (N), for coulomb 

2 

interaction with = 1.45. Multipole modes are denoted by corresponding I values and n = 1 
level indicate the breathing mode. 
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